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Recent astrophysical and terrestrial experiments have motivated the proposal of a dark sector with 
GeV-scale gauge boson force carriers and new Higgs bosons. We present a search for a dark Higgs 
boson using 516 fb~^ of data collected with the BABAR detector. We do not observe a significant 
signal and we set 90% confidence level upper limits on the product of the Standard Model-dark 
sector mixing angle and the dark sector coupling constant. 



PACS numbers: 12.60.-i,14.80.Ec 

While the astrophysical evidence for dark matter is 
now overwhelming, its precise nature and origin remain 
elusive. Recent results from terrestrial and satellite ex- 
periments have motivated the proposal of a new, hidden 
gauge sector under which WIMP-like dark matter parti- 
cles are charged [iHl]. An Abelian gauge field, the dark 
photon A', couples this dark sector to Standard Model 
(SM) particles through its kinetic mixing with the SM 
hypercharge fields [J]. In this framework dark matter 
particles can annihilate into pairs of dark photons, which 
subsequently decay to SM particles. The dark photon 
mass is constrained to be at most a few GeV to be com- 
patible with astrophysical constraints 0, |^ . In a mini- 
mal model , the dark photon mass is generated via the 
Higgs mechanism, adding a dark Higgs boson h' to the 
theory. The mass hierarchy between these two particles 
is not constrained, and the dark Higgs boson could be 
light as well. 

A consequence of this scenario is the possibility of 
probing a light dark sector at low-energy e^e~ collid- 



ers 0, Q and fixed-target experiments 0, [l^ . Searches 
for dark photon production have yielded negative results, 
and constraints have been derived on the mixing strength 
between the SM and the dark sector, e, as a function of 
the dark photon mass Q. 

The Higgs-strahlung process, e+e^ — s> A'h', h' — > 
A'A', might offer another gateway to a dark sector. This 
reaction is of particular interest, since it is one of the few 
process suppressed by a single factor of e, and the back- 
ground is expected to be small. If observed, this reaction 
could provide an unambiguous signature of physics be- 
yond the Standard Model. The event topology depends 
on the dark Higgs and dark photon masses. While Higgs 
bosons heavier than two dark photons decay promptly, 
their lifetime becomes large enough to escape undetected 
for nifi' < mA'- Moreover, the dark photon width is 
proportional to mA'f^, and its decay can be prompt or 
displaced, depending on the value of these parameters. 
At BaBAR energies, the decay length in the detector is 
O(IOO) /im or less for ruA' > 250 MeV and e > 10''', and 
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dark photon decays can be considered as prompt in this 
regime. 

We report a search for dark Higgs production in the 
Higgs-strahlung process. The measurement is performed 
in the range 0.8 < nih' < 10.0 GeV and 0.25 < ttia' < 
3.0 GeV with the constraint to^' > 2mA'- To avoid any 
experimental bias, the data are not examined before the 
selection procedure is finalized. The data sample used in 
this analysis consists of 521 fb^^ of data collected mostly 
at the T{4:S) resonance, but also including luminosity 
at the T{3S) and T{2S) peaks, as well as off-resonance 
data. A sample corresponding to ^ 10% of the data 
(optimization sample) is used to optimize the selection 
criteria and is discarded from the final dataset. This 
sample is treated entirely as background for optimization 
and background studies. 

The BaBAR detector is described in detail elsewhere 
[l3 | . Charged particle momenta are measured in a track- 
ing system formed by a five-layer double-sided silicon ver- 
tex detector and a 40-layer central drift chamber both 
immersed in a 1.5 T axial magnetic field. Electron and 
photon energies are measured in a CsI(Tl) electromag- 
netic calorimeter. Charged-particle identification (PID) 
is performed using an internally reflecting ring-imaging 
Cherenkov detector and the energy loss dE/dx measured 
by the silicon vertex detector and central drift chamber. 
Muons are mainly identified by the instrumented mag- 
netic flux return. 

Signal events are generated by MadGraph [l^ for 
about 40 different hypotheses of dark photon and Higgs 
boson masses. The hadronization of the A' qq {q = 
u, d, s, c) decay is performed by JETSET Jj^]. The detec- 
tor acceptance is studied using Monte Carlo (MC) simu- 
lation based on GEANT4 [l^ . Time-dependent detector 
inefficiencies, as monitored during data taking periods, 
are included in the simulation. 

The e+e~ A'h',h' A' A' reaction is either 
fully reconstructed in the 2{l'^l~)7r~^Tr~ and 

Z+Z^2(7r+7r^) final states {I ~ e,/i), or partially recon- 
structed in the 2(^+/i~) -I- X and fj.^iJ,~e^e~ + X chan- 
nels, where X denotes any final state other than a pair 
of pions or leptons. The 2(e+e^) + X mode suffers from 
significantly more background than the other channels 
and is excluded. The first modes are collectively referred 
to as "exclusive modes" , as opposed to "inclusive modes" 
for the 2{l^l~) + X channels. The inclusive modes are 
only considered in the region ttia' > 1.2 GeV, since their 
contribution is small below this threshold and the back- 
ground level becomes large. 

The event selection proceeds by first reconstructing 
dark photon candidates from pairs of oppositely-charged 
tracks identified as electrons, muons or pions by PID 
algorithms. In addition, the helicity angle of the 
electron in the dark photon rest frame, ag, must 
satisfy cosae < 0.9. The background from acciden- 
tal e+e~ pairs exhibits a peaking component near 



costte 1, while signal events are broadly distributed. 
Events are then processed according to the follow- 
ing sequence of hypotheses until a match is found: 
6^, 4^2e, 2/i4e, 6e, 4/i27r, 2^2e27r, 4e27r, 2/Lt47r, 2e47r, 4/z + 
X, 2/i2e -I- X. This order is chosen to minimize the 
cross-feed between channels and the efficiency loss due 
to misclassification. 

Additional criteria are applied to increase the purity 
of the signal. Exclusive modes must contain exactly six 
charged tracks, and the invariant mass of the three dark 
photon system must be larger than 95% of the e~^e~ 
center-of-mass energy. The dark photons are then fit- 
ted, constraining the tracks to originate from the in- 
teraction point. The fit probability is required to be 
larger than 10^^. Finally, the largest mass difference 
between the dark photon candidates, AM, must be less 
than 10 — 240 MeV, depending on the final state and the 
dark photon masses. The distribution of this variable 
after all other selection criteria are applied is displayed 
in Fig. [l]for the 2e47r final state. The signal peaks near 
AM ^ 0, while the background is concentrated towards 
higher values. 

Inclusive modes are selected by requiring two leptonic 
dark photon candidates with similar masses. The two 
dark photons are fitted, constraining the four leptons to 
originate from the interaction point. Events with a fit 
probability less than 10^^ are discarded. The remaining 
dark photon is then identified as the system recoiling 
against the two lepton pairs. The cosine of its polar angle 
in the laboratory frame must be less than 0.99 to remove 
radiative QED events. Finally, the masses of all dark 
photons must be compatible within their uncertainties. 

A total of six events are selected by these criteria: 
one 4/x27r, two 2/i47r, two 2e47r and one 4fi + X events. 
No candidate containing six leptons survives the selec- 
tion. The distribution of the dark photon mass versus 
the dark Higgs boson mass is shown in Fig. [21 Three 
entries, corresponding to the possible assignments of the 
decay h' — > A' A', are considered for each event. Be- 
sides the contribution of p — > vr+vr" or w ^ tt+tt^ 
decays near ttia' 0.7 — 0.8 GeV, no significant sig- 
nal is observed. This result is consistent with the two 
events observed in the optimization sample, assumed to 
be background. Given these limited statistics, a second 
background estimation based on the full dataset using 
same-sign combinations, such as (e"'"e~)(yLt+/i^)(/x~/i~) 
or (e+e+)(//~/i~)X, is used as a cross-check. Both meth- 
ods predict background levels consistent within their sta- 
tistical uncertainties. 

Using uniform priors in the cross-section, 90% confi- 
dence level (CL) Bayesian upper limits on the production 
cross-section are derived for each mode separately as a 
function of the dark Higgs and dark photon masses. The 
(m/i' , iTLA' ) plane is scanned in steps of 10 MeV in both di- 
rections between 0.8 < mh' < 10 GeV and 0.25 < tua' < 
3 GeV. For each mass hypothesis, the signal region is 
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FIG. 1: Distribution of AM, the largest mass difference be- 
tween the three dark photon candidates after all other selec- 
tion criteria are applied for the 2e47r final state. The data 
are shown for opposite-sign combinations from the optimiza- 
tion sample (plain circle) as well as an additional background 
estimation, described later, of same-sign combinations from 
the full dataset (open circle). The Monte Carlo prediction 
for m/j/ — 3.0 GeV and m^/ = 0.5 GeV is shown as a solid 
histogram. 



as signal every observed event in the signal region. The 
systematic uncertainties are included by convolving the 
likelihood of each final state with Gaussian distributions 
having variances equal to the systematic uncertainties 
described below taking correlations into account. 

The efficiency is determined for several values of dark 
photon and Higgs boson masses, and is linearly in- 
terpolated between the known points. The efficiency 
includes acceptance, trigger, selection criteria and the 
dark photon branching fraction. The branching frac- 
tions into leptons and hadrons are given by BF{A' ^■ 
£+£-) = 1/(2 + R), BF{A' hadrons) = R/{2 + R) 
and BF{A' tt+tt") = BF{A' hadrons)cr(e+e" 
7r+7r~)/cr(e+e^ — > hadrons), where R denotes the ratio 
^(e+e" hadrons)/cr(e+e^ — > ^i'^^i^) [11]. The effi- 
ciency increases from a few per mille in regions with small 
branching fractions to 33% for the six electron mode in 
the region mA' < 0.2 GeV. It drops rapidly in the region 
nih' < 0.8 GeV and mh' > 10 GeV, as tracks produced by 
dark photon decays have a low transverse momentum or 
are emitted close to the beam and are not reconstructed. 

The limits on each channel are then combined to ex- 
tract 90% CL upper limits on the e+e" — A'h',h' — 5- 
A' A' cross-section. The results are displayed in Fig. [31 
The limits are typically at the level of 10 — 100 ab. 
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FIG. 2: Distribution of dark photon mass (m^O versus the 
dark Higgs mass (m^/) for the final data sample. Three en- 
tries are plotted for each event, corresponding to the possible 
assignments of the decay h' — >■ A'A'. 




FIG. 3: Upper limit (90% CL) on the e+e" ^ A'h',h' 
A'A' cross-section as a function of the dark photon and dark 
Higgs masses. The limits in the lj- and i^-mesons regions are 
orders of magnitude larger than the average limits and are 
masked to avoid overflow. 



taken as the interval mh' — SfTm^, < < "^/i' + 3cr„i^, 
and vriA' - 5am^, < mA' < tua' + 3am^, , where cr™^, 
(cTm^, ) denotes the corresponding dark photon (Higgs) 
mass resolution. An asymmetric range is used to accom- 
modate the non-Gaussian tail of the low-mass side of the 
signal. The dark photon (Higgs) mass resolution varies 
between 2 — 17 MeV (3 — 55 MeV), depending on the dark 
photon (Higgs) mass and final state. While setting the 
limits we adopt the most conservative approach, treating 



The major contribution to the systematic uncertainty 
arises from the extrapolation procedure used to deter- 
mine the efficiency, which is estimated by comparing the 
extrapolated value to the nearest known point. This un- 
certainty increases from 1% to 8% in some corners of the 
phase-space. The uncertainty on the branching fractions 
ranges from a few per mille to 4%. The uncertainty due to 
the modeling of A' hadron decays in inclusive modes is 
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estimated by comparing different fragmentation models. 
This systematic is found to be 4% reflecting the hmited 
sensitivity of the selection procedure to the hadronic sys- 
tem produced by the dark photon decay. The uncertainty 
due to PID algorithms varies between 1.5% and 4.5%, 
assessed using high-purity samples of leptons and pious. 
Additional uncertainties include the determination of the 
track reconstruction efflciency (1.2%), luminosity (0.6%), 
and the limited Monte Carlo statistics (0.5% - 2.4%). 

The limits on the e+e~ A'h' , h' A' A' cross sec- 
tion are finally translated into 90% CL upper limits on 
the product aije^, where an = g\,/^iT and is the 
dark sector gauge coupling [7| . The results are displayed 
in Fig.|4]as a function of the dark photon (Higgs) mass for 
selected values of the dark Higgs boson (photon) mass. 
Values down to 10~^° — 10^^ are excluded for a large 
range of dark photon and dark Higgs masses. These re- 
sults assume prompt dark Higgs boson and dark photon 
decays. 

In conclusion, a search for dark Higgs boson production 
has been performed in the range 0.25 < ruA' < 3GeV 
and 0.8 < m/j/ < lOGeV for nih' > 2mA'- No signal 
has been observed and upper limits on the product of 
the mixing angle and the dark coupling constant in the 
case of a hidden sector with an Abclian Higgs boson have 
been set at the level of 10"^" — 10^^. Assuming ao — a, 
these measurements translate into limits on the mixing 
strength in the range 10^^ — 10^"^, an order of magnitude 
smaller than the current bounds. 
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FIG. 4: Upper limit (90% CL) on the product ao^ as a func- 
tion of the dark photon mass for selected values of dark Higgs 
boson masses (top) and as a function of the dark Higgs boson 
mass for selected values of dark photon masses (bottom). 
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